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Abstract A critical review of gravitational wave theory is made. It is pointed out that the 
usual linear approach to the gravitational wave theory is neither conceptually consistent nor 

■ mathematically justified. Relying upon that analysis it is then argued that — analogously 
! to a Yang-Mills propagating field, which must be nonlinear to carry its gauge charge — a 
I gravitational wave must necessarily be nonlinear to transport its own charge — that is, energy- 

^ • momentum. 
(D 

1 Introduction 

Although they have not yet been directly detected, there are compelling experimental evidences 
for the existence of gravitational waves, coming basically from the orbital period change of 
O^, binary pulsar systems [1]. These evidences, however, do not give any clue on their correct 

I form. This subject is actually plagued by theoretical difficulties that obscure their physical 

properties [2] . One of the main problems refers to the concept of energy localization [3] . On the 
one hand, the strong equivalence principle seems to forbid the energy of the gravitational field 
, to be localized On the other, from the point of view of field theory, it should be possible 

' to define an energy density for the gravitational field. Relying on this view, Synge says in the 

■ preface of his classic textbook [5] that in Einstein's theory, either there is a gravitational field or 
I there is none, according to as the Riemann tensor does not or does vanish. This is an absolute 
' property; it has nothing to do with any observer's world line. According to Synge, therefore, 

f-*) ■ the energy of the gravitational field should be localizable independently of any observer. In 

. the same line, Bondi argues that in relativity a non-localizable form of energy is inadmissible, 

' because any form of energy contributes to gravitation and so its location can in principle be 

found [6] . This controversial point is the origin of a long-standing discussion on the energy and 
^ . momentum transported by a gravitational wave. 

^ I The study of gravitational waves involves essentially the weak field approximation of Ein- 

stein equation, the field equation of general relativity. The geometrical setting of this theory 
is the tangent bundle0 at each point of spacetime, there is attached a Minkowski tan- 
gent space with coordinates x"" = This structure allows a clear distinction between 
coordinate system and reference frame [7j. Given a Lorentz frame 

e'^ = e%dx^ (1) 

it is possible to describe it in infinitely many coordinate systems {x^}. Observe that, although 
e" is coordinate independent, under a general spacetime coordinate transformation 

x'' ^ x'^" = x^" + (2) 



^We use the Greek alphabet {^,v,p, ... — 0, 1,2,3) to denote spacetime indices, and the first half of the 
Latin alphabet (a,fe, c, ... = 0,1,2,3) to denote algebraic indices related to the tangent Minkowski spaces, 
whose metric is assumed to be riab ~ diag (+1, —1, —1, —1). 
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its components e"^ change covariantly: 

(3) 

On the other hand, there are infinitely many frames e", differing each other by local Lorentz 
transformations, 

e'^'=A"'fee^ (4) 

which can be written in the very same coordinate system {x^}. Since a Lorentz transformation 
does not change the Minkowski metric r]ab, the spacetime metric 

S'At!/ — Vab 6 /^C jy, (5) 

is easily seen to be frame independent. As a consequence the Christoffel connection 

r^'^i. = \g''^ {d^gxu + dyQx^ - dxg^y) , (6) 
as well as the corresponding Riemann curvature tensor 

= - + t\^t\, - rvr^^, (7) 

are both frame independent. Although frame independent, the Riemann tensor components 
transform covariantly under the general spacetime coordinate transformation ([2]). 

2 Linear Approximation of Einstein Equation 
2.1 Wave Equation 

To study gravitational waves, one has first to obtain the weak field approximation of Einstein 
equation 

with Q^jj the source energy-momentum tensor. This is achieved by expanding the metric tensor 
according to 

9tiv = 'n^v + £ h(r)fiu + h(2)^u + • • • , (9) 

where e is a small parameter introduced to label the successive orders of this perturbation 
scheme. At this point, it is important to call the attention to the following fact. When the 
metric tensor is expanded according to ([9]) , we are automatically assuming that there is a back- 
ground Minkowskian structure in spacetime, with metric rj^^. Accordingly, the gravitational 
waves are interpreted as perturbations 

+ ... (10) 

propagating on that fixed Minkowskian background. This interpretation is consistent with 
general relativity, as well as with the point of view of field theory, according to which a field 
always propagates on a background spacetime [8j . 

Assuming expansion ([9]), the first order Ricci tensor is 

R{i)fj.v = dpT'^^-^^iy — di^T'^^^^p . (11) 
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Using the first order Christoffel connection 

^(i)fiu = \ri^^ {d^,h(^)Xu + duh(^)X^ - dxh^^)^^) , (12) 
the Ricci and scalar curvature tensors become, respectively, 

R(i)iiv = \ + d\d^h^^u + dxdyh^^)^ - d^d^h^^^^ (13) 

and 

= D/i,,) + dxd.h^,^^, (14) 

where □ = —rjP^dpdx is the flat spacetime d'Alembertian, and h^i^ = /i^^ja- In consequence, 
the first order sourceless gravitational field equation becomes 

° (^(i)Mi^ - Vfiiy hii)) - d^diyhfi) - r]f,y dad/sh'^f^ + d^dxh^y^y + dydxh^^^f, = 0. (15) 

Now, as is well known, wave equation (jlSp is invariant under general spacetime coordinate 
transformations. Analogously to the electromagnetic wave equation, which is invariant under 
gauge transformations, the ambiguity of the gravitational wave equation can be removed by 
choosing a particular class of coordinate system — or gauge, as it is usually called. The most 
convenient choice is the class of harmonic coordinate systems, which at first order is fixed by 

r?'^'^rf,)^, = 0. (16) 

In terms of the metric perturbation, this condition gives 

dphP^)f, = \dp\^y (17) 

In this case, the field equation (|15|) reduces to the relativistic wave equation 

□ = 0, (18) 

where we have introduced the new variable 

<i). = /iri)^-i'5^/i(i). (19) 

The harmonic coordinate condition ()17p can then be expressed in the Lorentz-like form 

d^r^,)^ = 0. (20) 

2.2 Linear Waves 

A plane-wave solution to the relativistic wave equation (jlSh has the form 

4)(i)^y = ^(1)/,!. exp[z/cpx''], (21) 
where = A(i)u^ is the polarization tensor, and the wave vector satisfies 

kp kP = 0. (22) 
The harmonic coordinate condition (j20p . on the other hand, implies 

^/x^fi)i^ = and (9^/ifi)iy = (23) 
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Still in analogy with the Lorentz gauge in electromagnetism, it is possible to further spe- 
cialize the harmonic class of coordinates to a particular coordinate system. Once this is done, 
the coordinate system becomes completely specified, and the components Af^^^i, turn out to 
represent only physical degrees of freedom. A quite convenient choice is the so called transverse- 
traceless coordinate system (or gauge), in which [9] 

Afi)p = and Af,^,Ul',^ = 0, (24) 

with an arbitrary, constant four-velocity. Since the trace condition A^^^p = implies that 

(pWfiu = h(i)p,u, we also have 

/ifi)P = 0, /ifi), = and A;^/ifi), = 0. (25) 

Now, although the coordinate system {x^} has already been completely specified (the 
transverse-traceless coordinate system), we still have the freedom to choose different local 
Lorentz frames e°. Since the metric Qf^u = rjab e-°'^e!'u is invariant under changes of frames, the 
metric perturbation will also be invariant. In particular, it is always possible to choose a frame 
fixed at one of the particles, called proper frame, in which C/^pj = . In this frame, as can be 
seen from the second of the Eqs. (j25p . 

/ifDO = (26) 

for all /i. Linear waves satisfying these conditions are usually assumed to represent a plane 
gravitational wave in the transverse-traceless gauge, propagating in the vacuum with the speed 
of light. Its physical significance, however, can only be determined by analyzing the energy 
and momentum it transports. 



2.3 Energy and Momentum Transported by Linear Waves 

The energy-momentum tensor of any matter (or source) field ip is always proportional to the 
functional derivative of the corresponding Lagrangian with respect to the spacetime metric. 
Since a derivative with respect to the metric does not change the order of the Lagrangian 
in the matter field both the Lagrangian and the energy-momentum tensor will be of the 
same order in the field variables ^. For example, since Maxwell's Lagrangian is quadratic in 
the electromagnetic field, the corresponding energy-momentum tensor will also be quadratic. 
Now, as is well known, the gravitational field is itself source of gravitation. This means that 
the gravitational energy-momentum density must appear explicitly in the gravitational field 
equation. In fact, the correct form of the wave equation (llSp is 

levrG 

° ^wtJ-i^ — ^wi^^ = 0, (27) 

with the first-order gravitational energy-momentum pseudotensor. As a consequence, 

in the linear approximation, the gravitational energy-momentum density is also restricted to 
be linear. Since the energy-momentum density is at least quadratic in the field variable, 
it vanishes in the linear approximation, yielding the wave equation (1180 . Of course, if this 
energy-momentum density is used to calculate the energy and momentum transported by 
linear gravitational waves, the result will obviously be that they do not carry neither energy 
nor momentum. 
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The consistency of this result can be verified by analyzing the generation of linear waves. 
In the presence of a source, the first-order field equation reads 



with Q{X)\iv the first-order source energy-momentum tensor. One solution is the retarded 
potential 

(x,t) = ^J e(,)^,(x',tO, (29) 

where the gravitational source is considered in the retarded time t' = t — \x — x'\/c. Now, as a 
consequence of the coordinate condition ([^0]) . it is easy to see that 

d.Qw^," = 0. (30) 

Instead of the usual covariant derivative, is conserved with an ordinary derivative at the 

first order. Since this is a true conservation law, in the sense that it leads to a time conserved 
charge — which in this case represents the source energy and momentum — we can conclude 
that in the linear approximation a mechanical system cannot lose energy and momentum in 
the form of gravitational radiation. This is consistent with the fact that linear gravitational 
waves do not transport energy and momentum. 

It is important to observe that, in the case of the electromagnetic field, the linearity of the 
wave equation does not restrict the electromagnetic energy-momentum tensor to be linear. To 
understand this point, let us write Maxwell equation in the form 

d^F^'''-f = r, (31) 

where we have introduced the electromagnetic self-current j'^ . As is well known, the linearity 
of electromagnetism restricts the electromagnetic self-current to be linear, and consequently to 
vanish: j'^ = 0. In fact, the electromagnetic wave is well known not to transport its own charge 
— that is, electric charge. The consistency of this result can be verified by observing that, 
in the electromagnetic case, the conservation law corresponding to (j3U|) is the conservation of 
the electric four current, 5^ = 0. Such conservation law precludes the electromagnetic wave 
to transport electric charge, but not energy and momentum. This is a fundamental difference 
between electromagnetic and linear gravitational waves. 

3 A Glimpse on Nonlinear Waves 

At the second order of an iterated perturbation scheme [H] , the harmonic coordinate condition 
reads 

7?^'^rf,,^, = 0. (32) 

This is equivalent to 

'■'ij.y^ (2)1' 

where 



S^'/'fa)^ = 0, (33) 



^(2)1/ — ^(2)1' ~ ^^'^'^ ^(2)) (34) 
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with /i(2) = hf(^2)p- III these coordinates, the second order gravitational field equation can be 
written in the form 

^ IGvrG ^ IGvrG ^ 

'-' '?(2)y t{2)y — ^{2)!^) (35) 

where t^^)^ = *{2)i^(^(i); ^(i)) represents all terms coming from the left-hand side of Einstein 
equation, in addition to the d'Alembertian term. It can be interpreted as the energy-momentum 
pseudotensor of the gravitational field. 

Now, as can be seen from Eqs. ()33p and (j35p . the second-order total energy-momentum 
tensor is conserved: 

^2). + %,u\ = 0. (36) 

The source energy-momentum tensor, on the other hand, as determined by the second order 
Bianchi identity, is conserved in the covariant sense: 

At the second order, therefore, the source energy-momentum tensor is not conserved. In fact, 
the above covariant conservation law is not a true conservation law, but simply an identity 
governing the exchange of energy and momentum between gravitation and matter [10]. As a 
consequence, in contrast to what happens at the first order, at the second order a mechanical 
system can lose energy in the form of gravitational waves. 

Far away from the sources, the second order gravitational waves are governed by the source- 
less version of the wave equation ()35p . which can be written in the form 

□ (^f^), = iV^(/i(i),/i(i)), (38) 

with A^/^j.(/i(i),/i(i)) = (167rG/c'^)tf2) jy. Since /i(i)^;/ is known, the problem reduces to that of 
solving a wave equation whose source is given. 

3.1 Emission of Gravitational Waves 

Considering the wave zone, that is, distances much larger than the dimensions of the source, 
the power emitted per unit solid angle in a given directiorH x = x/r, with r = \x\, is given 

by in 

^ = .= *'<*-°>. (39) 

where < t"^^ > is the energy flux averaged over a spacetime dimension large compared with 
1/uj, with uj the frequency of the wave. Now, in a perturbation scheme, different orders of 
the gravitational field expansion will give rise to different-order contributions to the power 
emitted. This means that we have to expand 

P = eP(i) +e2p,^, + ... . (40) 
At the first order, t'fX vanishes, and the corresponding power emitted will also vanish: 



= r^x'<t%>=0. (41) 



Middle letters fc, . . . = 1, 2, 3) of the Latin alphabet will be used to denote space indices. 
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This is consistent with the conservation law (I30p . and also with the fact that first-order (that 
is, linear) gravitational waves do not transport energy and momentum. 
At the second order, the power emitted by a mechanical system will be 

^ = (42) 

As already discussed, the energy-momentum pseudotensor depends only on the first-order 
solution. This means that Eq. (I42|) coincides with the usual expression found in all textbooks. 
It yields, therefore, the well known expression for the power emitted by a mechanical source, 
and in particular the usual quadrupole radiation formula. There is a difference, though: as 
a second-order effect, the energy and momentum emitted cannot be transported away by a 
linear wave, as is usually assumed, but only by a nonlinear gravitational wave. 



4 Effects on Free Particles 

4.1 The Geodesic Deviation Equation 

Let us consider, as usual, two nearby particles separated by the four- vector In the context 
of general relativity, this vector obeys the geodesic deviation equation 

Vc/Vc/r = iiV/3^''t^'e^ (43) 

where = dx^/ds, with ds = g^^ dx^dx'^ , is the four-velocity of the particles. Now, each 
order of the gravitational field expansion 

R°'iiui3 = £ R'(i)iiuf3 + R%)tivi3 + ■ • • , (44) 

which follows naturally from ([9]), will give rise to a different contribution to For consistence 
reasons, therefore, this vector must also be expanded: 

r = Co)+eeri)+£'C2) + -- - • (45) 

In this expansion, represents the initial, that is, undisturbed separation between the par- 
ticles. Of course, as the four-velocity depends on the gravitational field, it must also be 
expanded. We then write 

U^ = Ui^,,+eUi^,,+e^U^,, + ... , (46) 

where C/jq) is a constant arbitrary four-velocity, which depends on the choice of the initial con- 
dition — or equivalently, on the choice of the local Lorentz frame from which the phenomenon 
will be observed and measured. Using then the freedom to choose this frame (see section YI7^ . 
we can choose a frame fixed at one of the particles — called proper frame. At the lowest order, 
therefore, the four- velocity ?7(q) can be expressed in terms of the observer proper time S(o), that 
is, [/(Q) = dx^/ds^o)j where 

dsfo) = V^^•y dxf'dx'' (47) 

is the fiat spacetime quadratic interval. Since, in the proper reference frame, the particles are 
initially at rest, we have that 

C/(^)^ 5^0 = (1,0,0,0), (48) 
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which means that, in this frame, the proper time S(o) coincides with the coordinate 

Before proceeding further, it is crucial to observe that the geodesic deviation equation (j43p is 
covariant under general coordinate transformations. In consequence, the equations emerging at 
each order of the perturbation scheme will also be covariant up to the corresponding order. Due 
to this covariance, one can either solve the equation in a general coordinate system and then 
specialize to the transverse-traceless coordinates, which are the coordinates used to describe 
the gravitational wave, or specialize first to these coordinates and then solve the ensuing 
differential equations. Here, we will adopt the first alternative, which we consider far more 
convenient. 

4.2 Linear Effects 

At the lowest order, the geodesic deviation equation (H3]) becomes 



"^(0) 



(0) 



0. (49) 



This is an expected result, in the sense that ^[^^ is simply the undisturbed separation between 
the particles. Without loss of generality, we can take the solution to be ^[^^ = constant. At 
the next order, the geodesic deviation equation becomes 



""'(0) 



+ C/(o) dp (r"j/3-y C/(o)) ^[^) — R'^i^fj.u/B ^^(0) U^o) ^(0)- (50) 



Substituting U^y^ as given by Eq. (fiS]) . we get 



ds 



+ c^or")/3o Cm - R(i)00i3 S,^o)- (51) 

Using then the first-order Riemann tensor 

-R"i)/ii//3 = duT^^^^p - OpV^^^fj,^, (52) 

it reduces to 

Canceling SoF^j^q C(o) on both sides, we get 

§^ = -«ooefo), (54) 

UOf. 



'(0) 

where 



rfijoo = I rf' (2 ao/i(i)po - ap/i,i)oo) . (55) 

Specializing now to the transverse-traceless coordinate system, where the components /i(i)po 
vanish identically, we obtain 

« ?(i) 



0. (56) 
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This means that, in the hnear approximation, the particles are not affected by gravitational 
waves. 

For the sake of completeness, let us consider also the other usual approach to the same 
question, which consists in calculating the proper distance between two nearby particles. Such 
a distance is defined by 

Al = J \ds^\'/\ (57) 

with ds'^ = g^udx'^dx'^. Due to the fact that the gravitational wave solution is written in 
a very specific coordinate system, this computation has to be done either in this specific 
coordinate system or using a coordinate-independent formulation. Here, we will follow the 
second alternative^] Using then the metric expansion ([9]), and keeping terms up to first order, 
we obtain 

ds^ = ds^Q^ + e dx'^dx'^ , (58) 

with ds^Qj given by Eq. (j47p . Choosing again the proper reference frame fixed at one of the 
particles, we can write 

ds^ = {l + e [/(^) U^,;) ds% , (59) 
where = dx^/ds^oy However, in the transverse-traceless gauge. 

Up to first order, therefore, the proper distance between two nearby particles is simply 

A/ = y" 1^41^2, (60) 

which means that it is not changed by a first-order gravitational wave. Only at the next order, 
where nonlinear effects will be in action, will the proper distance be affected by a gravitational 
wave. 

Although some of the components of the first-order curvature tensor R^^^^up do not vanish 
in the transverse-traceless coordinate system, as we have seen, they are unable to induce 
any movement on free particles. Since these components are different in different coordinate 
systems, it is difficult to get some insight on their meaning. A better way to understand the 
effects of curvature is to inspect the invariants ^12] constructed out of R'([^^uf3i as for example 
the Kretschmann 

and the pseudo-scalar 

t' — , p"/'/^'^ p — i,Q/3Ap p IJ-^ n 

with -k denoting the Hodge dual. As an easy calculation shows, both invariants vanish every- 
where for the linear solution (j2ip . It is not surprising, therefore, that such waves are unable 
to impart energy and momentum to free particles. 

^In all approaches found in the literature, this computation is made in a coordinate-dependent formulation. 
However, these approaches are not consistent, in the sense that different coordinate systems are used to describe 
the wave and the position of the particles. In fact, whereas the gravitational wave is described in transverse- 
traceless coordinates, the position of the particles are usually specified in Cartesian coordinates. 
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4.3 Some Comments on the Usual Procedure 



In all approaches found in the literature, the effect produced by linear gravitational waves on 
free particles is computed by first specializing the geodesic deviation equation to the transverse- 
traceless coordinates, and then by solving the corresponding differential equations. Of course, 
since the coordinate system is fixed from the very beginning, these approaches are not mani- 
festly covariant. Moreover, all of them present some unclear points, which we pass to discuss. 

It is usually claimed that, by choosing a frame fixed at one of the particles, the Christoffel 
connection can be made to vanish, not only at one point, but along the whole world-line of a 
comoving observer fixed at this framed In other words, it is usually claimed that [3] 

= aor-,,, = 0, (61) 

where use has been made of the fact that, in this frame, the proper time S(o) coincides with the 
coordinate x^. This result is then used to eliminate the connection-term from the left-hand 
side of Eq. (|53|) . but not the same connection-term appearing in the right-hand side of the equa- 
tion [14J. The argument used to justify this procedure is that, since R^^^^up is a tensor, it can 
be calculated in any coordinate system, not necessarily in the transverse-traceless coordinates. 
This argument is, however, misleading. The reason is that, once the coordinate system has 
been completely specified, the equations become not manifestly covariant, and consequently 
any argument based on covariance cannot be applied. Actually, consistency reasons do require 
that, since the geodesic deviation equation is written in the transverse-traceless coordinates, 
the Riemann tensor components on the right-hand side of this equation must necessarily he 
expressed in the same coordinate system. As a consequence, if the connection-term in the 
left-hand side of Eq. (|53|) is found to vanish, the connection-term appearing in the right-hand 
side of the same equall^n must also be assumed to vanishl Furthermore, observe that, since 
the left-hand side as a whole is also a tensor, the use of the same criterion should preclude the 
elimination of its connection term. 

The important point, we repeat, is that arguments based on covariance cannot be used 
when the coordinate system has been completely determined. Notice, for example, that in the 
transverse-traceless coordinate system, we obtain from Eq. (j52p the relation 



Rl.^up Uj^,, C/ro) = K) 9p (rri)/37 U;,) ■ (62) 

Of course, this identity is not inconsistent with covariance because it is valid only in that 
specific coordinate system. Therefore, when the geodesic deviation equation ([5T]) is restricted 
to the transverse-traceless coordinates, the above two terms can be consistently canceled from 
the geodesic deviation equation ([5T]) , yielding the equation of motion ([56]) . 



^ We notice in passing that there is an apparent problem with this claim: as we have discussed in section [T] 
the Christoffel coefficients do not depend on the frame, and hence cannot be made to vanish by a choice of a 
frame. In order to make a connection to vanish along a curve, it is necessary to transform to a locally inertial 
coordinate system, sometimes called Galilean coordinates [T^. However, this is not possible here because the 
coordinate system has already been completely specified. Remember that the gravitational waves are described 
in the transverse-traceless coordinate system. In addition, as a simple inspection shows, not all components of 
the Christoffel connection (|12p vanish in the transverse-traceless coordinates, which means that this coordinate 
system is not Galilean. 

^Observe that this would not mean that all components of the Riemann tensor vanish: only its projection 
along the particle four-velocity would vanish, that is, -R(\),j,y/3 U^g^ (7^^) = -R(\)00/3 = 0. As can be easily verified, 
other components of the Riemann tensor would still be non-vanishing to comply with its covariant behavior. 
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5 Final Remarks 



Considering that the sources of gravitational waves are at enormous distances from Earth, 
the amphtude of these waves when reaching a detector on Earth are expected to be so smah 
that the hnearized theory is usually assumed to be enough to accurately describe them. In 
other words, to comply with the idea that the radiative solutions to the gravitational field 
equation should satisfy a linear wave-equation, these waves are usually assumed to carry not 
enough energy and momentum to affect their own propagationl^l However, these assumptions 
are clearly misleading. The basic reason is that this is not a matter of approximation, but a 
conceptual question. Either a gravitational wave does or does not carry energy: if it carries, 
it cannot satisfy a linear equation. It is, therefore, conceptually unacceptable to assume that 
a gravitational wave satisfying a linear equation is able to transport energy and momentum. 
Only a nonlinear gravitational wave is able to do it. The above assumption would correspond 
to assume that, if a solitary- wave solution to the Korteweg-de Vries equation had a small- 
enough amplitude, its evolution could be accurately described by a linear equation. Needless 
to say, this is plainly wrong. 

The problem of defining an energy-momentum density for the gravitational field has a close 
analogy with the problem of defining a gauge current for the Yang- Mills field, whose sourceless 
field equation is |15| 

Q^pA„^ + j^-" = 0. (63) 

The piece 

represents the gauge pseudocurrent, which is present only for non-Abelian gauge groups (for 
which f^BC 7^ 0). Due to the explicit presence of the connection fj_, this current is clearly 
not gauge covariant, in analogy with the gravitational energy-momentum pseudotensor, which 
is not covariant under general coordinate transformations. A fundamental characteristic of 
the Yang-Mills field is that it carries its own charge, as for example the color charge carried 
by gluons in chromo dynamics. This simple property precludes the existence of Yang-Mills 
linear waves. In fact, linear Yang-Mills waves — that is, solutions of the linearized Yang-Mills 
equation — are waves that do not carry their own charge, a property also revealed by the 
fact that, in the linear approximation, the pseudocurrent j^'^ vanishes. The Yang-Mills field, 
therefore, is essentially nonlinear — otherwise it is not a Yang-Mills field. 

Let us consider now the gravitational field, whose sourceless field equation, in the so called 
potential form, is written as [16] 

d>.iV^Sxn - ^ V^tx' = 0, (65) 

where g = det{gf^^) and Sx^f^ = —S\^^p is the superpotential. This equation presents a structure 
similar to the Yang- Mills equation ()63p . with the gravitational energy- momentum pseudo- 
density {8ttG /c'^) htx^ playing the role of self-current. This similarity comes from the fact 
that, analogously to Yang- Mills fields, the gravitational field also carries its own charge — 
in the case, energy and momentum. Because the energy-momentum pseudotensor is at least 
quadratic in the field variables [8], it vanishes in the linear approximation. Up to first order, 

^This is a commonplace in gravitational wave theory. For a textbook reference, see [11] , 
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therefore, the sourceless field equation ()65p becomes the field equation (llSp . which in harmonic 
coordinates reduces to 

□ = 0. (66) 

In this form, it is similar to Maxwell's equation in the Lorentz gauge. This could be interpreted 
as an indication that, in this approximation, the gravitational field loses its similarity to Yang- 
Mills fields, becoming similar to the electromagnetic field, whose self-current also vanishes. 
However, for the same reason a Yang-Mills field must be nonlinear to carry its own charge, a 
gravitational wave must be nonlinear to transport energy and momentum. Otherwise, it is not 
a gravitational wave. It should be remarked that, besides providing the dispersion relation, 
the solution /i(i)^,/ of the linear equation ()66p is physically meaningful, in the sense that the 
higher-order solutions depend on it. Nevertheless, it does not represent by itself a gravitational 
wave which, like a Yang-Mills propagating field, is a strictly nonlinear phenomenonlll 

Summing up, although the use of a perturbative analysis could suggest the existence of 
linear gravitational waves, physical reasons show that a gravitational wave is essentially a 
nonlinear phenomenon. Of course, any mechanical system — like a binary pulsar, for example 
— can lose energy in the form of gravitational radiation. The usual expression for the power 
emitted by a mechanical source, and in particular the quadrupole radiation formula, give a 
correct account of this energy. What our results say is simply that this is a second-order effect, 
and as such the energy and momentum cannot be transported away by a linear wave, but by 
a solution of a nonlinear wave-equation. At the lowest order, it is given by the second order 
wave equation 

□ ./.f,), = iV^(/i(i),/i(i)). 

The traveling wave solution to this equation, therefore, represent the physical gravitational 
wave in the sense that it is able to transport energy and momentum, and consequently to be 
detected. It is the wave to be searched when looking for gravitational waves. An analysis of 
these waves, therefore, as well as of their effects on a test particles, turns out to be crucial for 
determining their signature on gravitational wave detectors. 
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